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Abstract: Astrophysical objects such as active-galactic nuclei (AGN) and gamma-ray 
bursts (GRBs) can be sources of high energy, astrophysical neutrinos. The decay of charged 
pions produces electron and muon-flavor neutrinos from the primary decay of the pion, and 
from the secondary decay of the resulting charged lepton. At low energies we expect the 
flavor ratio <& Ve : <& u : <& Vt to be 1 : 2 : at the source. We are interested in the flavor 
ratios as measured on Earth after the neutrinos propagate over cosmic distance scales from 
the source. If we only consider vacuum flavor transition probabilities between the three 
active flavors then we expect a measured flavor ratio of 1 : 1 : 1 up to small corrections 
from #13 and non-maximal #23- When we include mixing with two additional flavors of 
sterile neutrinos then we see corrections to this ratio up to ~ 30%. Furthermore, if we 
consider energy-loss of the charged leptons involved in the pion decay from cosmic sources 
then these flux ratios depend on the neutrino energy. We examine these energy-dependent 
flavor ratios using a specific model for sterile neutrino mixing, and compare to expected 
ratios when only three neutrino flavors are considered. At energies E v > 1 TeV the flavor 
ratios observed in experiments such as IceCube can probe the existence of sterile neutrinos. 
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1 Introduction 

High-energy, charged pions are produced by p7 or p-nucleon scattering in astrophysical 
objects. The charged pions may be from gamma-ray bursts (GRBs) [1-3] where they are 
produced with a flux oc E~ 2 [3], and also they may result from gamma ray-jets in 
active galactic nuclei (AGN) [2, 4, 5]. These pions produce neutrinos through the process 
7r + — > v „ + fj> + — > Vn + v e + e + + Vn, so we expect a 1 : 2 : flavor ratio from the source 
where all ratios are given as : <& v : & Ut . Neutrino propagation is affected by flavor 
oscillations, and oscillation probabilities are included in the calculation of final flavor states. 
Therefore measured neutrino fluxes on the Earth are flavor- mixing model dependent. When 
we only consider three flavors of neutrinos, we expect a flavor ratio 1 : 1 : 1 as measured 
on the Earth [6, 7], up to corrections from a non-zero #13, and non-maximal #23 [8]. 

We must consider cooling effects from astrophysical sources on the charged leptons and 
pions which produce these neutrinos. Energy losses factor in from astrophysical sources 
such as losses due to synchrotron radiation [1, 3] or inverse-Compton emission [3], as 
in GRBs, or adiabatic cooling as in AGN [3, 5]. These losses give the measured ratios 
a neutrino-energy dependence. Different models of energy-losses lead to different flavor 
ratio profiles; measurements of these neutrino flavor ratios over several decades of neutrino 
energies may allow one to determine the source of the neutrinos depending on the profile. 

In this paper we will consider the expected neutrino flavor ratios when couplings to 
two additional flavors of sterile neutrinos are included. Corrections to flavor ratios when 
one additional flavor of sterile neutrinos is included have been considered before in [7, 9]. 
In section 2 we will discuss different models for the cooling effects of the charged pions 
and leptons. In section 3 we will describe the formalism of neutrino flavor-oscillations, and 
briefly describe the model of sterile neutrino mixing used here. Finally, we will compare 
the predictions from a five-flavor case to the three-flavor case of the Standard Model in 
section 4. 
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2 Models for charged pion/lepton energy loss: source effects 



We now want to consider different models of charged pion production as well as energy- loss 
for charged pions and charged leptons prior to their decay. 

We consider an astrophysical source which produces charged pions with a flux & w oc 
E~ k . Energy losses are also assumed to follow the power law E x = dE x /dt oc —E^. If 
we consider pions produced by the process p + J — > n + tt in a GRB then we have a 
flux of pions which behaves as $> n oc E~ 2 [1, 3]. The presence of strong intergalactic 
electromagnetic fields lead to radiative energy losses of the pion and subsequent charged 
leptons prior to decay. Synchrotron energy loss, and likewise inverse-Compton emission, 
correspond to n = 2 [3]. If we consider adiabatic energy loss as a result of expansions of 
the 7r + / fi + plasma then this corresponds to n = 1 [3]. In the process p + 7 — > n + tt, 
although the neutron can decay to produce P e , the fraction of the proton's energy that this 
neutrino carries is much smaller than the fraction carried by the neutrinos resulting from 
pion decay. We will therefore ignore the contribution to fluxes from the resulting neutron. 

The energy dependent fluxes of v^, v e and z7„ at the source are calculated for example 
in [3], and the results are given here. For the production of f M from the primary decay of 
7r + we have 

*i(E„) = S „(-s„) (1 - fc)/n e- s "r ( — , 0, -s n ) (2.1) 
For the production of and v e from the secondary decay of the /i + we have a source flux 

= ^-r(-^) (1 - fc+n)/n {<T*T (^,0, -sn) 

q I \ n J ( 2 _ 2 ) 

V-V^r(^,o,- 9 " Sn )} 

where s n = \ ^ ^4^°°' \ > Q = ^ M ' c °°| , and r(^p,0,— s n ) is a lower-incomplete gamma 
function. E XyCOO \ is the energy at which the time, T^cooi to achieve significant energy loss 
(cooling) due to synchrotron radiation or adiabatic energy loss is the same as the decay 
time, r^decay [3]. Energy loss from astrophysical sources changes quickly with energy, such 
that r^cooi oc E x /E x , and so r^cooi/T^decay « E~ n , and also (E WtCOO i/E^ coo i) ~ 10 2/n [3]. 

In the next section we will use the energy-dependent neutrino flux at the source to 
calculate fluxes on the Earth of the three active neutrinos when couplings to two sterile 
neutrinos in the 3+2 MM are considered. 

3 Neutrinos oscillations/propagation effects and the 3+2 minimal model 

Neutrino flavor-states are linear combinations of the different mass eigenstates. Given N 
flavors of neutrinos, there are N neutrino mass eigenstates, and mixing is determined by a 
unitary N x N matrix, U. Generically this mixing is given by 

| fa) = ^ U aiWi) 
i 
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3+2 MM 


lAm^K eV 2 ) 


lAm^KeV 2 ) 


\UeA 


Pa] 


l^4l 


IM 


045 


NH 


0.47 


0.87 


0.149 


0.127 


0.112 


0.127 


1.87T 


IH 


0.9 


1.61 


0.139 


0.122 


0.138 


0.107 


1.47T 



Table 1. Results of the fits for the 3+2 MM for both mass hierarchies. 



where U is parameterized by several mixing angles and phases; given TV neutrino flavors 
there are N(N — l)/2 mixing angles, and up toN(N + l)/2 phases. In the case of the 3v of 
the Standard Model the matrix U is given by the PMNS matrix, parameterized by current 
global fits given in [10]. 

We can compute oscillation probabilities for the process v a — > v$ ; probabilities depend 
on the neutrino energy as well as the length of propagation, L, and are given by 

P a p(L/E) = 5 aP - 4 MKPpiUpjUcd) sin 2 

i>j 

+2J2lm(U* i U^ i U* pj U aj )sin 

i>j 

We are neglecting potential terms which arise from neutral-current and charged-current 
weak interactions with matter, and only considering vacuum oscillations. Matter potentials 
are proportional to the density of electrons or nucleons in the propagating medium, and 
these densities are very low in most astrophysical cases which makes matter potentials 
negligible. Some astrophysical objects, such as the fireballs from GRBs are sufficiently 
larger that only very light neutrinos are affected, therefore it is a valid choice to exclude 
any matter interactions. For example, in a GRB fireball the electron density is ~ 10 10 — 
10 12 cm -3 , and so the MSW resonance would only significantly affect neutrinos with mass 
differences less than Am 2 ~ 10~ 12 eV 2 [11]. 

Since we are considering neutrinos from astrophysical sources we have a very large 
propagation length, and therefore we are in the limit x = Am^L/AE S> 1. In this limit 
oscillations are very rapid, and oscillation terms in the probability take on their average 
values, sin 2 (x) — > \ and sin(2:c) — > 0. Therefore the probabilities in the flux calculations 
become 

P a p = (PaftiL/E)) = 5 a p -2Y / Re(U* i U /3i U^U aj ) = ^ \U ai \ 2 \U Pl \ 2 (3.2) 

i>j i 

The model for sterile neutrino mixing used here involves two sterile neutrinos, with 
two additional sterile mass eigenstates where 7714.5 ~ C(eV 2 ) [12]. This model is a minimal 
extension of the Standard Model because it involves adding only two Standard Model gauge 
singlet Weyl fields; each field corresponds to a unique sterile flavor. The additional terms 
in the Lagrangian are given by 

Cbsm = -l a L Y a ^u J R - ^4 C M|4 + h.c. 

where index summation is implied, and c indicates charge conjugation. Y is the 3x2 
Yukawa matrix, <3? is a scalar, SSB field, and Mr is a diagonal, 2x2 mass matrix. A 



AE 



Am^L 
2E 
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basis is chosen for this model such that the first mass eigenstate becomes massless, and 
Am| 2 = Am§!. 

The details of the parameterization and fits for this model can be found in the original 
paper [12]. To summarize, the 3+2 MM is parameterized by four mixing angles (the 
three angles of the Standard Model and one additional angle which mixes the two sterile 
mass eigenstates), three phases which includes a relative phase between the sterile mass 
eigenstates, and four non-zero mass eigenstates. The results of the parameterization are 
summarized in table 1. The authors in [12] found that the normal hierarchy provides a 
better fit to neutrino flux data which includes anomalies, and therefore we will work in the 
normal hierarchy (NH) for the remainder of this paper. 

4 Active neutrino flavor ratios on the Earth 

As the neutrinos are allowed to propagate from the source to the detector on Earth there 
will be transitions between flavor states. The final states measured on the Earth are given 
by 

P 

where transition probabilities used here are given by the limiting case in eqn. 3.2. We will 
be examining the case where the pions at the source are produced in a gamma-ray burst, 
and the flux goes as Q n oc E~ 2 (k=2). Plots of the flavor ratios for different models of 
pionic, and leptonic energy losses can be seen in Fig. 1(a), and the comparison to the three 
flavor case of the Standard Model can be seen in Fig. 1(b). 




Figure 1. (a) A plot of $„ Q ,/<& I , e . The solid curves consider fluxes where the ir + and /i + energy 
losses are due to synchrotron radiation or reverse-Compton emission (n=2); the solid red curve 
(upper) is a = fi, while the solid blue curve (lower) is a = t. The dashed curves consider fluxes 
where the tt + and /i + energy losses are adiabatic (n=l); the upper dashed curve is a = fi, while 
the lower dashed curve is a — r. Flux from tt + and tt~ decay are considered here, (b) The ratio 
($y Q /$yJ 3+2 MM /($„ Q /<£>„„ ) 3!y . The red (upper) is a — /j,, and the blue (lower) curve is a = r. 
The ratios are calculated for n=2 only. 

We can see the flavor transition in the neighborhood around E v = E^ coo \ where the 
flux ratios transition from ~ 1 : 0.9 : 0.7 to ~ 1 : 1.4 : 1.1. Determining E^ coo \ is model- 
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dependent. An interesting feature can be seen in Fig. 1(b), where the ratio &uu/^v e 
predicted by the 3+2 MM is larger than the ratio predicted by transitions between the 
three flavors of the Standard Model over all decades of neutrino energies, but especially at 
high energies \og w (E v /E^ CQO {) > 0.5. This implies larger probabilities of P(y a — > v^) via 
some process such as v a — > v s — > involving one or more transitions from active to sterile 
back to active flavors of neutrinos. This predicts that we would expect the measured 
flux to be larger at ultra- high energies, and the ratio Q Vli /§ Ue is significantly higher with 
two sterile neutrinos than the same quantity for only three flavors. This would provide a 
clear test for the presence of sterile neutrinos in the 3+2 MM scheme. Flavor ratios as a 
function of neutrino energy when E^ iCOO \ = 1 TeV can be seen in Fig. 2(a); £^ iCO oi = 1 TeV 
approximately corresponds to a model of gamma-ray bursts associated with the collapse 
of a massive star where cooling is a result of inverse-Compton emission [3]. 




500 1000 1500 2000 500 1000 1500 2000 

E v [GeV] Ey [GeV] 

(a) (b) 

Figure 2. (a) A plot of the same quantity as in Fig. 1(a), where Sp.cooi = 1 TeV. This muon 
cooling energy approximately corresponds to a model of gamma-ray bursts associated with the 
collapse of a massive star [3]. (b) Comparing the ratio of tracks (iO to cascades (v T + v e ) between 
the 3+2 MM and the prediction with only three flavors where E^_ coo \ = 1 TeV. 

For energies E v < 1 TeV it is not possible to separate electron neutrinos from tau neu- 
trinos. In detectors such as IceCube [13] muon neutrinos produce tracks whereas electron 
and tau neutrinos contribute to cascades at these low energies. At high energies it may 
be possible to differentiate between electron and tau flavors by observing lollipops or u T 
charged-current double-bangs [14]. The comparison of the tracks to cascades ratio between 
the 3+2 MM and three flavors can be seen in Fig. 2(b), where the energy region prohibits 
the differentiation between the tau and electron flavors. 

5 Conclusions 

Examining neutrino flavor ratios as a function of the neutrino energy can potentially allow 
one to determine the source of the neutrinos. Shapes corresponding to k = 2 and n = 2 as 
in the solid curves of Fig. 1(a) may indicate a GRB source, whereas shapes corresponding 
to k = 2 and n = 1 as in the dashed curves may indicate an AGN source. Calculating flux 
ratios also has the potential for determining whether there are flavors of sterile neutrinos. 
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The flux predictions in this paper also provide a method for determining the existence 
of sterile neutrinos, specifically whether we have a minimal extension of the Standard Model 
involving two sterile flavors. If GRBs are associated with the collapses of massive stars 
then £?ucooi < 1 TeV due to inverse-Compton emission [3], and therefore experiments such 
as IceCube [13] which can measure high-energy neutrinos over several energy decades are 
good candidates for measuring these flavor ratios. 

Verifying these flavor ratio profiles requires statistics from many neutrino events. 
Nearby gamma-ray bursts which provide more than 1-2 neutrinos are infrequent (~ 
1/century), and GRBs typically have durations of 0.1 - 100 seconds [2] so the number of 
neutrino events from a single GRB will be small; it is necessary to take measurements from 
many GRB sources. Although AGN emit high energy radiation over a period lasting weeks 
[2] the number of neutrino events expected from a single AGN is low. In both GRB and 
AGN cases we must wait a long time in order to measure enough events to extract flavor 
ratios from background. 
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